Light localization signatures in backscattering from periodic disordered media 
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The backscattering line shape is analyticaUy predicted for thick disordered medium films where, 
remarkably, the medium configuration is periodic along the direction perpendicular to the incident 
light. A blunt triangular peak is found to emerge on the sharp top. The phenomenon roots in the 
coexistence of quasi-lD localization and 2D extended states. 

PACS numbers: 42.25.Dd,42.25.Hz 



Introduction. Since the mid-eighties [ij the coherent 
backscattering (CBS) has been one of the pilots of stud- 
ies of Anderson locaUzation of light. Indeed, manifesta- 
tions of weak localization (WL) in the CBS line shape 
have been well documented for various disordered dielec- 
tric media 0, while how strong localization (SL) affects 
CBS has been a long term fascinating subject [Sj]. The 
last decade has witnessed spectacular progress on CBS 
near 0, Q or far below Q the localization transition, 
which undoubtedly is an intellectual challenge both ex- 
perimentally and theoretically. Indeed, to prepare strong 
scattering media and to extract localization from medium 
absorption are highly restrictive Q, while the failure of 
perturbation theory Q-crucially mapping the pictorial 
reciprocal paths into (diagrammatical) one-loop approxi- 
mation -enforces the invention of a nonperturbative 
theory to allow a microscopic analysis. Despite of these 
difficulties a common belief is that SL is responsible for 
rounding the CBS sharp top 0,11,11,11]. 

In studies of light localization much attention has been 
paid to fully disordered media. There have also been in- 
creasing interests on other medium structure such as dis- 
ordered photonic crystals [sl where the Bloch symmetry 
is slightly destroyed by impurities, and systems with per- 
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feet Bloch symmetry [9| where WL is analytically found. 
Most interestingly, the recent invention of so-called pla- 
nar random laser introduces a novel medium structure 
[lo| . which consists of a random gain layer sandwiched 
by two mirrors. It was then conjectured that SL in the 
layer plane might be responsible for the laser emission 
[lO|. To prove it is yet a nontrivial task which may be 
traced back to the striking feature of the partially dis- 
ordered structure. Indeed, (from the geometrical optics 
view) two perfect reflection mirrors map the medium to 
an extended periodic one but fully disordered inside the 
primitive cell (apart from the mirror symmetry ). This 
immediately raises many important questions. For exam- 
ples, does Anderson localization exist in such structure? 
If so, can it be probed by CBS measurements ? 

Unfortunately, the interplay between Anderson local- 
ization and the periodicity is a difficult issue. In general 
the common wisdom regarding localization may be dras- 
tically modified and very little has been known, among 
which are: a new scale essential to WL [1| may ap- 
pear, and whether the constructive interference between 
reciprocal paths encompasses WL depends on periodic 
medium configuration [HI . 

In this letter these two problems will be investigated 
for a simplified model-thick disordered medium film with 
a dielectric function periodic in the x-direction (Fig. [T]). 
The lattice constant is a and the primitive cell consists of 
randomly positioned point-like scatterers filling the half 
space z > . The film section is uniformly illuminated by 
a beam of stationary unpolarized light (with the wave- 
length A and the frequency 17) perpendicular to it. Mi- 
croscopic analysis is provided for the angular resolution 
of light intensity near the inverse incident direction for 
sufficiently large times. 

Qualitatively, as shown in (Fig. [2]), the traditional line 
shape [1, 2] develops at \q±\ > lixja and is sharpened at 



1 < 



\q\\ < 27r/a, eventually a blunt triangular peak 



FIG. 1: Light backscattering from a periodic thick disordered 
medium film. Inset: the film section. 



emerges at \q^\ <C ■ Here q± = (27r/A)sin0 and 
^0 = navD with ly the photon density of states at f2 
and D = 1/3 the bare diffusion constant {I the transport 
mean free path and the velocity c ~ I). Quantitatively, 
analytical predictions are made for A ^ Z ^ a <C ^ e'/'*' , 
where the last inequality ensures photon states to be far 
from 2D SL In particular, the line shape a{6) is 



2 



singular at \q± \ = , 27r/a , around which 

r (i + /|g^|)-2^/o(0), k±|>^; 

{l + l{0)}lo{0), 0<^-\q^\«^; (1) 
I {1 + l/Co)- a/ (271^0) l\q±\, \q±\ ^ 

while a smooth interpolation between the last two lines 
is expected. Here the enhancement factor 'y{d) = 
(31/2^0) [I -a\q^\/{2Tr)]. 

Qualitative picture. The well known Bloch theorem 
allows us to reduce the photon motion into an effec- 
tive one within a primitive cell dictated by the Bloch 
wave number kp . The enhanced backscattering finds 
its origin in the constructive interference-described by 
low-energy hydrodynamic modes-between two counter- 
propagating photons (so-called cooperon) each of which 
carries a Bloch wave number kp ,k'p , respectively, kp + k'p 
plays the role of "Aharonov-Bohm flux" . The gauge in- 
variance then requires that the transverse (x-direction) 
hydrodynamic wave number 2'kN / a with N an integer 
satisfies q±^ = 2nN/a — {kp + k'p) . Hydrodynamic modes 
with TV ^ (iV = 0) describe 2D (quasi-lD) motion. 

As kp,k'p G [— 7r/a,7r/a) there are two contributions 
which correspond to two successive N responsible for the 
line shape. For > 27r/a both iV ^ . Therefore, 
these two hydrodynamic modes are inhomogeneous in the 
transverse direction, and extended in the longitudinal (z- 
) direction because of a ^ / e'/"^ . As usual the diffusive 
2D motion then leads to the traditional line shape. 

For < 27r/a the hne shape is contributed by both 
2D diffusive and quasi- ID motion since the two succes- 
sive integers are now 1 (or —1) and 0. The former (or 
latter) occupies a portion of a\q±\/2Tr (or 1 — a\q±\/2Tr). 
In the quasi-lD geometry there are two scales: I and • 
An incident flux decays over the scale ~ I then diffuses 
inside the medium and eventually exit at z ^ I . The 
penetration length is ^ \q±\^^ ^ I . 

For \q±\ closed to 27r/a photons penetrate into the 
medium of a distance <C via quasi-lD diffusive motion. 
Upon penetration they may self-intersect then propagate 
around the formed loop along the same direction-so- 
called diffuson (Fig. [3] (a)). The probability is larger in 
quasi-lD than in 2D. As a result two initially closed but 
counter-propagating photons have a larger probability to 
be brought back to their starting point and the backscat- 
tered light intensity is thereby enhanced. At z ~ / the 
probability of forming such a loop is ~ l/D. Conse- 
quently, the inverse diffusivity at z ~ Z increases with an 
amount of l/{avD) ^ l/£,a ■ Noticing that the backscat- 
tered light intensity is proportional to the inverse dif- 
fusivity [l| and taking into account the weight of 
quasi- ID motion, we find that the traditional line shape 
is magnified by a factor oi 1 + ^{0) . 

Notice that the quasi-lD cooperon and diffuson have 
different masses: D{kp + k'pY — Dq\^ and D{kp - k'pY , 
respectively. As \q±\ decreases the diffuson tends to ac- 
quire a larger massive and be damped. Consequently, all 




FIG. 2: The predicted (solid) versus traditional (dotted) line 
shape (symmetric with respect to = Q). The parameters 
are a = 5Z, = 10/ (setting / = 1). Inset: the blunt trian- 
gular peak at \q±_ \ <C ~ 0-1 ■ 

the constructive interference involving diffuson-cooperon 
couphng (e.g.. Fig. [3] (a)) tends to be suppressed. Op- 
posed to this the cooperon becomes less massive. Conse- 
quently higher order loop-wise interference paths involv- 
ing solely cooperon-cooperon couphng (e.g.. Fig. [3] (b)) 
accumulate and eventually dominate the backscattered 
light intensity at ^ ■Co'^ ; where photons penetrate 
deeply into the medium forming SL in the bulk. 

In this region from one-parameter scaling hypothesis 
we expect that the diffusion coefficient exponentially de- 
cays from the interface, i.e., — lnZ3(z) oc z/^o- There- 
fore, the average inverse diffusivity of the boundary layer 
increases also by an amount of ~ Z/^o • Thus for jg^l ^ 
^q"^ the quasi- ID motion contributes to the line shape 
(1 -t- //^o)[l — a|9_L|/(27r)] , together with the portion con- 
tributed by the 2D extended motion: /o(6')[a|(7j_|/(27r)] « 
a\q±_\/ {2t:) , and leads to a blunt triangular peak. 

General formalism. We then outline the proof [l^ . Let 
us start from the retarded (advanced) Green's function: 
G^i^(R,R') = (R|{V2-f [l-|-e(R)]±iO+}-i|R') de- 
scribing the propagation of the electric field. The fiuc- 
tuating dielectric field e(R) has zero mean and vanishes 
for z < , while for z > is periodic in x and satisfies 
(e(R)e(R')> =^5{v- v') 5{z - z') Y.Nez 5{p - p' - 
Na) , where (3m e Z) p = a; — ma G [0, a) stands for the 
relative x-coordinate in the primitive cell. The scattering 
is elastic with I = An/ A . Using the effective medium ap- 
proximation the incident field is i?(R) = EQe^^^'^~''/^' 
[l3|. Alternatively, one may introduce the Green's func- 
tions: g^i^ir,r';kp) = (rK^i^ _ H{kp) ± iO+)-^\r') 
for the effective motion within the primitive cell with 
Hikp) = - [{dp + ikpf + dl + dl] ~ 6(r) (r s (p, v,z)). 
The two sets of Green's functions are related through 

G«i^(R, R') = e^'-^^^-^'^g^iv, r'; kp) . (2) 
The albedo characterizing the radiation intensity in the 
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direction s = (sin 6* , , — cos 9) [l3| generally depends on 
the time t . With the single scattering event ignored the 
albedo at t ^ oo , denoted as a{9) , can be shown to be 



dRidR2 e ' i^n^ (R-i i -R-2 ) (R-2 j R-i 



g^as.(R.-RO (r^^ r^) (Ri, R2)] (3) 



in the unit of A^Sq/IGvt with (■ • • ) the average over the 
fluctuating dielectric field e(R) . The first (second) term 
gives the background intensity ao (line shape a{9)). No- 
tice that the y-dependence of Green's functions is irrel- 
evant and will be ignored from now on. Substituting 
Eq. ([2]) into Eq. ([3]) gives (up to an irrelevant overall 
normalization factor) 



ao 



JJ^ dzdz'e-'^Y.yoiz^z'-kp.kp), (4) 



a{9) = U dzdz'e-'^Y.Y."'yNiz,z';kp,k'^),{5) 

'0 AT U Ut 



where we have quantified the propagators: diffuson and 
cooperon introduced above to be >'°(r, r'; fcp, /c^) = 

gg,(r,r';fcp)g,-^,(r',r;fcp and y^ir^-kp^k';) = 
tjj?j^2 (r, r'; fcp) ^/^2 (r, r'; fcp , and introduced their 
Fourier transformations, i.e., 3^^'^(r, r'; fcp, fc^ = 
E^ez e*(P-p')2-A^/a yY{z, z'; fcp, fc^) . The partial 

summation Y.k„k', = '^k,M'^ 4p + fe; ,2^Ar/a-g^ • 

One may proceed to sum up over all the (maximally 
crossing) ladder diagrams (e.g., Refs. 0, [i3|) producing 
a bare diffuson (cooperon) 3^'^(^) . Nevertheless such 
a diagrammatic expansion fails in the nonperturbative 
analysis-intrinsic to localization. Instead, for interac- 
tionless systems such as photons to fulfill this task the 
supersymmetric method-to which we switch below-turns 
out to be perfectly suitable T]. Conceptually, the intro- 
duced slow varying Q-field interprets the bare diffuson 
and cooperon as Goldstone modes of the spontaneous 
supersymmetry breaking, and encapsulates their mutual 
interactions underlying localization through the nonlin- 
ear constraint: = 1 . Technically, it can be shown 



y"^^(r,r';fcp,fc;; 



27 



{stT[kA+T-Q{r)A-T^kQ{r')]) (6) 



following standard derivations [3], where (5(r) — 
T~^(r)AT(r) is a matrix field with full orthogonal sym- 
metry, i.e., T G [/(2,2/4)/C/(2/2) x U{2/2). The av- 
erage is defined as (. . . ) = J D[Q]{- ■ ■)e~^^'^^ with 
Q\z=-zo = A, (5|p=o = Q\p=a , where the action F[Q] = 
^ ^^dz /;dpstr[(9,g)2 + [DpQf] with Dp = dp + 

i[fcT3, •] . str is the supertrace and all the matrices follow 
the definitions of Ref. • In addition, A± = (1 ± A)/2 , 
= (1 ± T3)/2, and the matrix k = diag(/cp , fc' ) is 



diagonal in the retarded-advanced sector. The bound- 
ary condition Q(— zq) = A accounts for the fact that the 
albedo is contributed by optical paths not crossing the 
trapping plane located at z = — zq ~ — 0.7Z « 
which coincides with the medium boundary. 

Eqs. (H])-® constitute the general formalism of 
calculating the albedo. By making an appropriate 
global rotation for the Q-field reflecting the symmetry: 
y^{v,v']kp,kp) — y'^{r,r';kp,~'k'p) , it can be shown 
that, similar to the fully disordered medium with con- 
served reciprocity the background intensity ao = 
a(0) . We then turn to analyze the line shape. 

Line shape. Three regions: (i) \q±\ > 27r/a, (ii) 
^0^ ~ |9J-I < 27r/a and (iii) \q^^\ < will be stud- 
ied separately. It can be shown that in (i) and (ii) 
Q{r) mildly fluctuates around A. A perturbation the- 
ory near it generates leading diffusive motion-described 
by the bare diffuson (cooperon)-and loop-wise interfer- 
ence (Fig.[3|) encompassing WL. Opposed to this, in (iii) 
though locked at A at the boundary Q{r) may strongly 
fluctuate along the z-direction driving photons into SL 
states, while remains homogeneous in the p-direction . 
Consequently the perturbation theory breaks down. 

In Eq. ([5]) the summation over N picks up two terms 
with successive integers N = N' , N' + 1 . In (i) they both 
do not vanish corresponding to 2D low-energy extended 
motion. The leading order perturbation then gives 



~d,Dd, + D 



2nN 



+ (kp T k'p) 



yN 'kpi k ) 



^■6{z-z'), y^'^'u^o^o- (7) 

Substituting its solution into Eq. ^ recovers the first 
line of Eq. (P) [Ii, [H. In (ii) and (iii) |A^| = , 1 . Tak- 
ing into account the weight of 2D motion, i.e., a|q'_L|/(27r) 
the term with |7V| = 1 is found to be 2^ 1^(9) . The 
term with = arises from the quasi-lD motion suf- 
fering from localization effects, and below will be calcu- 
lated separately for (ii) and (iii) with the simplified ac- 
tion F[Q] = 2^/o°°dzstr{(5^Q)2 + [iIts^Q]^} (since 
Q does not depend on p). 

For (ii) one-loop expansion show that Eq. ^ (N — 0) 
still holds reflecting the flux conservation. However, 
the diffusion cocfhcient acquires the position-dependence 
and is found to be D{z) = D{1 - (^d/cp ± Vl)"^ [1 - 
g-2z|fep±fcp/ ^j^jj ^j^g _j_ gjgj^ corresponding to y^ and 

yQ , respectively. Therefore, we justify microscopi- 
cally the crucial conjecture-position-dependent diffusion 
coefficient-made in Ref. Q at the one-loop level. Nev- 
ertheless the difference should be stressed that at the 
medium boundary z — the bare diffusion constant is 
protected against one-loop WL. Such an important prop- 
erty persists up to higher order loop corrections enforcing 



D(z = 0) = £». 



(8) 
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(a) 


















FIG. 3: One-loop (a) and typical two-loop (b) interference 
picture underlying quasi-lD weak localization. 



details are to be reported elsewhere llTl.) and (phe- 
nomenological) hydrodynamic methods [l8|. Observing 
that aie) « ^'EivEfc„fc/^^^(^^;fcp>^p) ^e need to 
consider only the returning probability-like propagator 
y^il,l;kp,k'p) in (iii). 

Applying the heat kernel method 16] to Eq. ([6]) we 
succeed to calculate 3^q"(z, z; 0, 0) exactly provided that 
T is lowered down to the GMat(3, 2|A) symmetry jT^]. 
The most important feature of the solution is the ex- 
ponential divergence 17[: (All the numerical factors are 
unimportant and not given here.) 



(10) 



In the bulk: z ^ \kp ± kp\^ the diffusion coefficients 
become homogeneous: D{z) — > D[l — (^ol^p ±^p|)"^] but 
strongly depend on their infrared cutoffs. By contrast, 
since Eq. ([5]) suggests that all the interfering optical paths 
contributing to the line shape reside in a boundary layer 
of size ~ \q±\~^ (> , for < 27r/a - \q^\ < ir/a the 
local diffusion coefficient of is simplified as D{z) w 
D{1 — 2z/^o) , and well approximated by D{z) « D/{1 + 
2z/^o) as the leading order Z/^o correction concerned. 
Replacing D in Eq. ^ with the latter we obtain: 



1 



[/_ (a;> )/+ (a;< ) - C /_ (x> )/_ (a;< )] , (9) 



where 2:>(<) = Cok_L|[l + ^ "i^x (min) {z , z'}] , C — 
f+{^o\q±\)/f-iCo\q±\) and f+{x) = xli{x) and f-{x) = 
xKi{x) with Ii{x) ,Ki{x) the modified Bessel functions. 
Eq. (O is then inserted into Eq. (O. For S,o\q± \ » 1 one 
may use the asymptotic expressions of the Bessel func- 
tions and eventually find the line shape given by the sec- 
ond fine of Eq. ([!]) [15|. At \q±\ ~ 7r/a interfering optical 
paths significantly penetrate into the bulk and suffer from 
stronger WL, resulting in a larger enhancement factor. 

For (ii) higher order expansion of Eq. ([5]) shows that 
the line shape, indeed, is contributed by interfering opti- 
cal paths forming the loop-wise structure classified into: 
diffuson-cooperon coupling where two paths may trace 
some loops in the same direction (e.g, Fig. [3] (a)), and 
cooperon-cooperon coupling where two paths trace all 
the loops in the opposite direction (e.g. Fig. [3] (b)). As 
shown below it is the latter leading to SL in the bulk at 
|9±| ~ ^cT^ : and the line shape at \q±\ ^ £,q^ is mainly 
responsible for by (radiative) SL states with lower sym- 
metric T . 

For (iii) to calculate the lineshape namely 
y^{z,z';kp,kp) exactly is a very hard task. Due 
to the broken translational symmetry the technique 
of Ref. breaks down and the solution there is not 
applicable. Instead, our method combines the mi- 
croscopic yet advanced mathematical theory namely 
performing super- Fourier analysis [l^ for Eqs. © and 
([S]) (This is far beyond the scope of this letter and the 



for z ^0 • On the other hand, because it was shown 
that bulk SL states display hydrodynamic behavior 
we expect in the presence of vacuum-medium interface 
the exponential divergence Eq. pO)) to be reflected at 
the same macroscopic level. To achieve this we notice 
that the boundary leakage introduces the level broaden- 
ing r(z) scaling as — Inr(z) z/^o On the physical 
ground the contribution of local currents to the restoring 
force-leading to SL-exponcntially decays in time, modi- 
fying the VoUhardt-Wolfle model [18] to be 



— Siz-z'), 



(11) 



with J^Q^I^^o = and the covariant derivative V = 
{dz , dp + iq±) , where the overall coefficient of the restor- 
ing force term is fixed by Eq. ([8]). Solving Eq. pT|) indeed 
confirms Eq. ^TU\i {qj_ = and the localization length 

e = eo). 

The presence of nonvanishing kp k, —k'^ alters the mi- 
croscopic symmetry and therefore the localization class. 
Indeed, (for |g^| < ^q"^) one may follow Ref. [ll|], av- 
erage Eq. ([6]) over \kp — k'p\ and subsequently obtain 
an effective action: F[Q] = dzstT{{d^Q)^ + 

[iq±A/2,Q]'^} . The T-field symmetry is lowered down 
to {/(I, 1/2)/;7(1/1) X (7(1/1) giving ^ = 2^0 0. 

On the other hand, the common belief of one- 
parameter scaling hypothesis (e.g. Ref. 20]) implies that 
the microscopic symmetry enters only through the lo- 
calization length leaving the hydrodynamic model unaf- 
fected. To find y^iz, z'; kp, k'p) wc insert the z-dependcnt 

level broadening: r(z) = £1/(2^0)^6"^ [13] and C = 2^o 
into Eq. (fTT|) . For t ^^/D the steady distribution is 
approached solving: —D{dze^^ dz — \q±\^e~^}yQ — 
{ttv)^^ 6{z — z') , which coincides with the diffusive model 
for the single channel SL [5]. Substituting the solution 
and the bare propagator, Eq. ([7]) into Eq. ([5]) one may 
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find 



m = 1- 



(1- 



3-v^l + (4eo9±)2 



a/(27r^o) l\q± \ , 



27r 

-1 
• 



(12) 



The first line liolds for \q^ \ ^ £,q^ and, apart from tlie fac- 
tor: [1— a|g^|/ {2tt)] accounting for tlie weiglit of quasi-lD 
motion, the first term resembles the rounded line shape 
of fully disordered media below localization transition [5| . 
Remarkably, the second line suggests that the line shape 
displays a blunt triangular peak in a very narrow region: 
\q±\<t:^o^ (inset of Fig. E]). 

Conclusions. Analytical studies of the coherent 
backscattering line shape have been presented for peri- 
odic thick disordered medium films which arrest both ex- 
tended and quasi-lD localization states. The result is ex- 
pected to be qualitatively correct for A < Z < a (< le''^'^) 
which, together with the realization of the perfect peri- 
odicity along one direction, may be well within the reach 
of up-to-date experimental conditions [ilOj . However, to 
study realistic media the present theory still needs to be 
extended so that the large size (namely the film thickness 
much larger than the lattice constant) effects and the par- 
ity (mirror) symmetry are taken into account, which is 
left for future work. 

I am grateful to S. Hikami for useful conversations, 
especially A. Altland, L. Zhou and M. R. Zirnbauer for 
important discussions. Work supported by Transregio 
SFB 12 of the Deutsche Forschungsgemeinschaft. 
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